Generation of twin Fock states via transition from a two-component Mott insulator to 

a superfluid 
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We propose the dynamical creation of twin Fock states, which exhibit Heisenberg limited inter- 
ferometric phase sensitivities, in an optical lattice. In our scheme a two-component Mott insulator 
with two bosonic atoms per lattice site is melted into a superfluid. This process transforms local 
correlations between hyperfine states of atom pairs into multi-particle correlations extending over 
the whole system. The melting time does not scale with the system size which makes our scheme 
experimentally feasible. 

PACS numbers: 03.75.Lm, 03.75.Mn, 03.75.Dg 
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Degenerate atomic Bose- and Fermi-gases [lfl provide 
an excellent starting point for engineering almost-pure 
many-particle quantum states which are an essential re- 
source for novel quantum technologies. This is illustrated 
by the realization of a Mott insulator (MI) in the low- 
est Bloch band of an optical lattice which can serve as a 
quantum memory However, new methods are neces- 
sary for attaining some of the most important TV-particle 
states, in particular those which allow a sensitivity en- 
hancement from the standard quantum limit oc 1/y/N 
to the Heisenberg limit oc 1/JV in quantum metrology 
Q. Here we propose a method for engineering twin 
Fock states [H starting from a two component MI with 
two bosonic atoms per site in an optical lattice d, 
Our scheme manipulates the hyperfine states a and b 
of atom pairs pinned to single lattice sites and decou- 
pled from one another in the MI regime. In this limit 
the dynamics is reduced to a set of identical two par- 
ticle problems and correlations between states a and b 
are accurately induced using Raman laser pulses or mi- 
crowaves and collisional interactions controlled via Fesh- 
bach resonances 0] ■ These local two-particle correlations 
are then transformed into multi- particle correlations ex- 
tending over the whole system by melting the MI into 
a superfluid (SF). This quantum melting can be exper- 
imentally implemented by adiabatically ramping down 
the depth of the lattice potential Q as shown in Fig. [Ha). 
By an appropriate choice of correlations created in the MI 
phase a twin Fock state emerges in the resulting two com- 
ponent superfluid. Using this state in a Mach-Zchnder 
interferometer (MZI), shown in Fig. [ijb), one can a P~ 
proach sensitivities scaling as 1/N. 

Our method is motivated by theoretical and experi- 
mental evidence 0, 0] that long range correlations build 
up quickly when melting a MI to a SF. Almost adia- 
batic melting can be achieved with ramping times t r 
on the order of a few tens of N/MJo where Jo is the 
typical tunneling amplitude between lattice sites during 
the melting process and M is the number of lattice sites 
Q. Importantly, in our scheme M oc N thus preventing 
t r from scaling with N. The achievement of sensitivi- 




FIG. 1: (a) Melting of a two-component two-atom MI into 
a SF. (b) The internal states a and b represent the arms of 
an interferometer with rotations Rl, R2 and phase-shift cj> 
induced by a Ramsey pulse sequence; <f> is measured by M. 



ties at the Heisenberg limit relies acutely on the relative 
number of atoms entering each port of the interferome- 
ter. This stringent requirement rules out simple schemes 
involving either a 7r/2 rotation into a two-component 
Bose-Einstein condensate (BEC) or a rapid splitting of a 
single-component BEC in a double-well potential, since 
both result in a binomial relative number distribution. 
While number squeezing to a twin Fock state could be 
achieved in the latter case by an adiabatic splitting, this 
timescale t r oc TV since M = 2 is fixed In contrast, 

the smaller t r for our scheme is less demanding on the 
suppression of dissipation which is known to degrade the 
achievable sensitivities 13, II | • We note another promis- 
ing route to generating a twin Fock state through the 
coherent dissociation of a molecular BEC as suggested in 

El. 



Our starting point is a MI with two atoms in each lat- 
tice site i given by | ^ a b) — IL I a ^)i @> E3- This state, 
with exactly the same number of a and b atoms, can 
be created by collisional interactions involving an auxil- 
iary state c as experimentally demonstrated in We 
analyze the melting of | ^ a b), and also of the superposi- 
tion state \ ^aa+bb) = Y[ii\ aa )i + | 66> i )/v / 2 obtainable 
from | ^ a b) by applying a 7r/2-Raman pulse. For this we 
use the two component Bose-Hubbard model describing 
the dynamics of atoms trapped in the lowest Bloch band 
of a sufficiently deep optical lattice. The corresponding 



2 



Hamiltonian is (h = 1) [2] 

(».i> « 



where Oi(6i) is the bosonic destruction operator for an 
a(&)-atom localized in lattice site i, hi — a\a,i and 
h\ = St Si, while denotes summation over nearest- 

neighbors. The parameter J a n,) is the tunneling matrix 
element for atoms in state a(b); V a tb) an d U are the on- 
site intra- and inter-species interaction matrix elements 
respectively. For simplicity we only consider the sym- 



metric case with J. 
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V. The ratio be- 



tween the matrix elements is determined by the lattice 
depth and additionally U, V can be controlled via 
Feshbach resonances 0| or by shifting the a and b atoms 
away from each other using state-dependent lattices 
We denote the total number of atoms in state a(b) as 



N. 



a(b) 



= and introduce Schwinger boson opera- 

nt? 



EMh + b]ai)/2, J y = EMbi ~ &}bi)/2i, 
-N h )/2 and J = (J x , J y , J z ). 



tors J x — 
X. = (N a 

We illustrate the outcome of melting two-component 
MI ground states with fixed total particle number N — 
2M assuming adiabatic evolution in an M-site sys- 
tem with periodic boundary conditions. The melting 
starts deep in the MI regime, where J can be ne- 
glected. For U < V the nondegenerate MI ground state 
is | ^ab) and this adiabatically melts into the nonde- 
generate SF ground state | * tf ) = \N/2) Aq \N/2) Bq oc 

(ioV^OBoV^kac) for V/J 0. Here all N/2 
a-atoms are in the same delocalized symmetric mode 
Aq oc Yli^ii a ^ N/2 b-atoms are in mode Bq oc J2i b\ 
and | vac) is the vacuum state. Thus, adiabatic melting 
of | ^ab) provides a direct means of obtaining a twin Fock 
state | \l/tf) with zero relative atom number difference. 

In the opposite case U > V the system exhibits 
spatial separation of the a and b components and has 
a large number of degenerate ground states. For an 
even number N a of a-atoms this degeneracy is lifted 
by completely connected hopping 14J and the state 

|*£p) = HU( l <N a /2)\ aa ) l Il(N a /2<j<N)\ bb )j} is adi- 
abatically connected to the nondegenerate SF ground 

state | N a ) A | N — N a ) Bg , where $ denotes symmetriza- 
tion over lattice site configurations. Our second ini- 
tial state | ^aa+bb) is a binomial superposition of ground 
states | ^^p) and therefore results in the melted state 
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After rotation, the overlap of this state with 
the macroscopic superposition state (IBj | ^ mn ) = 
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FIG. 2: The amplitudes |C m | 2 (left-bottom axes) and corre- 
sponding relative phase probability distribution P(0) (right- 
top axes) after the rotation Rl of (a) | ^tf) and (b) | ^ mac ) 
with N = 20. 



(\N) Ao I°)b + l°U \ N )b )/V2 is found to be = 
K*ma X |cxp(i7rJ a /2)|* mac )| = ^8/9 > 0.97 for N -» oo 
and this limit is monotonically attained with O > 0.97 
for N > 20. 

We can use either of the final states obtained by adi- 
abatic melting to realize sensitivities oc 1/N in the MZI 
setup shown in Fig. [TJb). For the twin Fock state 
| ^tf) Rl implements the conventional beam-splitter op- 
eration exp(i7r J x /2). In the case of | V^mac) Rl rotates 
the state according to exp(i7rJ a /2) as discussed above. 
Then a relative phase is induced in one of the arms 
of the MZI. The operation Rl can be implemented by 
rapid resonant 7r/2-Raman pulses, while the relative 
phase (f) could be induced by an appropriate off-resonant 
pulse. The achievable phase sensitivity inside the inter- 
ferometer (after Rl) can be computed from the relative 
phase probability distribution P{&). For a given state 

N 



I V>) = E m =o C m I to )a I N _ m ) B a tnis is computed as 
P S (A8) = \ EZ=o C ™eM-™iA0)\ 2 /( S + 1) where A8 
is a multiple of 2ir/(s + 1). The distribution P{9) is ob- 
tained by multiplying P S (A6) with (s + 1)/2tt and tak- 
ing the limit s — ► oo [11] . As shown in Fig. Ela) P(6) 
is sharply peaked for the rotated | ^/ t f) and Fig. [5{b) 
displays the oscillatory behavior of P(8) for the rotated 
I t mac ). The sensitivity 8<j) of each state is quantified by 
the half width at half maximum (HWHM) of P(9) around 
some fixed phase typically taken as 8 — 0. From this def- 
inition it can be shown for both states that 8(f> oc 1/N 
and so scale at the Heisenberg limit 

[lElGjl- The final 
step in the MZI consists of rotation i?2 and measure- 
ment M. of a phase dependent quantity which exhibits 
this Heisenberg limited phase resolution. For instance 
the measurement of parity as discussed in [lj| or J 2 as 
proposed in [ijj can achieve this. These methods exploit 
nonlinearities caused by atomic collisions and are thus 
realizable in our setup. 

We now analyze the achievable sensitivity for incom- 
plete melting at finite (residual) values of V/J before 
considering experimental imperfections arising from non- 
adiabatic ramping and particle loss. We restrict our con- 
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siderations to the measurement of j| and to the twin 
Fock state due to its greater experimental feasibility. We 
study these effects in terms of the noise Acf> on the cf> de- 
pendent observable J 2 . Error propagation theory gives 
A0 = AJl/\d{Jl)/dcj)\ where (J 2 ) and A J 2 are the av- 
erage and the spread respectively. For states which are 
zero eigenvectors of J z , such as | \&tf ) , this gives 



A0 : 



2 _ sin 2 0((j4) 



(J 2 ) 2 )+cos 2 <f>(J x J*J x ) 



4 cos 2 0(J 2 )' 



(2) 



-1/2 



At (j> = the sensitivity reduces to A(j>(0) = ^(J x , 
which can be expressed entirely in terms of the one- 
particle density matrices = (ajaj) and p\- = (b\bj) 

using (J 2 ) = \J2M + &) + y.Jl'ilt + t.c.)]/4. The 
initial MI state | ^f a b) with no off-diagonal correlations 
Pij = Pij = $ij yields the standard quantum limit Ac/) — 
1/V2N (cf. Fig. da) at V/J > 1). The final SF | * sf ) 
with long-range correlations p° = p \^ = 1 asymptoti- 
cally recovers the Heisenberg limit A<fi =1/ \/N 2 /2 + N 
(cf. Fig. |3fa) at V/J = 0). Thus the scaling Ac/) oc N~ a 
changes from a = 1/2 to a = 1 during the melting. 

The presence of a residual intra-species interaction V 
results in quantum depletion of the populations in the Aq 
and Bq modes which reduces the attainable sensitivity. 
We consider this effect for U = 0, achieved for example 
by fully separating the sites of a spin-dependent lattice. 
In the SF regime J>l^a translationally invariant sys- 
tem is well described by the A-conserving Bogoliubov 
wavefunction 



|^og>~(Aa) JV/2 (A6) Jv ^|vac) 



(3) 



where A a = (Aj 2 — S q ^o c q^q^-q) the q quasi- 

momentum modes for a, and similarly A& is defined in 
terms of Bq for component b. Within this ansatz the 
depletion is identical for both components and is given 



b y C = Eq^0"q = C q/( 1 _ C q) Where ^ = (4q^q) = 

(BqBq) and the amplitudes c q <C 1 can be solved in 
terms of V/J and are given in [l^|. The sensitivity 

is then A<j>(0) = l/[iV 2 (l - C) 2 /2 + N + 2 E q#0 K^- 
For a three dimensional (3D) optical lattice the deple- 
tion C = (V/J) 3 / 2 /3tt 2 v / 2 and the sum £ q #o n q = 
(Vyj)3/2( 37r _ 8)/l2V2 are constant in AT. Within the 
region of validity of the Bogoliubov ansatz there will thus 
be no degradation of the sensitivity scaling a. Similarly 
we find no decrease of a in 2D within the Bogoliubov 
ansatz. In contrast, for a ID system the depletion in- 
creases with since there is no true condensate for non- 
zero V and long-range correlations decay algebraically. 
Consequently the scaling a decreases with increasing V/J 
as shown in Fig. [U(a) for a finite number of particles. 

We explore the sensitivity scaling between the SF and 
MI limit in ID numerically using the Time-Evolving- 
Block-Decimation algorithm [19( . For specific values of N 




FIG. 3: (a) Scaling a in A(j) ~ N a against V/J of ground 
states of H with U — 0. The scaling was obtained from 
calculations with N = 20 up to N = 60 with M = N/2 for 
both the ID numerical results (o) in a box system and the 
iV-conserving Bogoliubov results (solid line), (b) Sensitivity 
A<j> as a function of V/J for N = 60. The dotted lines show 
exact results for the twin Fock SF (bottom) and the MI (top) 
states, while the thick solid line shows the Bogoliubov result, 
(c) A4> as a function of N for V/J = 20 (MI) and V/J = 1/2 
(SF) giving a of the corresponding points in (a). In all cases 
the dashed lines are to guide the eye. 



we compute the ground states of the Hamiltonian Eq. ([TJ) 
over a range of V/J and U — and calculate their sen- 
sitivities. The result for N — 60 in Fig. [^b) shows a 
smooth transition between the ideal SF and MI results 
and agrees with the Bogoliubov result for V/J < 2. By 
repeating this calculation for different N we extract the 
sensitivity scaling a. For all values of V / J the numerical 
results are well approximated by a power-law A(f> ~ N~ a 
as shown for the extreme cases in Fig. [3]Jc) . We combine 
these results in Fig. [3£a) which displays the scaling a as 
a function of V/J. The value of a decreases smoothly 
from the Heisenberg limit in the SF regime to the stan- 
dard quantum limit in the MI regime. As also shown 
in Fig. the lD-Bogoliubov ansatz Eq. is consis- 
tent with the numerics for V/J < 2 and the values of 
N used here. In the MI regime a number conserving 
particle-hole ansatz [2(| underestimates long-range cor- 
relations and therefore does not predict an increase of a 
from 1/2 at finite V/J. This emphasizes the importance 
of the growth of long range correlations for increasing a 
above the standard quantum limit. While our numerical 
calculations are performed for sizes of ID optical lattices 
currently used in experiments [3] we expect the scaling 
a to drop from a = 1 to a = 1/2 at vanishingly small 
V / J when N — > oo because of the lack of true long range 
correlations. In 2D and 3D this reduction is instead ex- 
pected to occur at finite V/ J close to the phase transition 
point (V/J) cr it below which long range correlations exist 
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FIG. 4: (a) Energy difference per atom (E — E )/NJ between 
| ^(tr)} and the final SF ground state | *Po} (left axis, dashed 
lines) and overlap \ {^(t r ) \ ^o) \ (right axis, solid line) as 
a function of the ramping time t r . Starting from the MI 
ground state at Vo/J = 20 we ramped V(t) down to a final 
Vt r /J = 1/2 using different ramping times t r and U = 0, 
N = 50 and M = 25. (b) Scaling a of Scj> with N for | * tf ) 
under the presence of atom-loss at a rate k obtained from 
results for N — 20 up to N = 60. Note that a does not reach 
unity because of the finite N used. 



at large N. 

In experiments the ramping time t r will be finite and 
so the melting is never perfectly adiabatic. We numer- 
ically compute the dynamical ramping of V/ J for the 
Hamiltonian Eq. fl} with U = from the MI ground 
state. In Fig. |4ja) the energy difference per atom be- 
tween the final ramped state \ ^(t r )) and the SF ground 
state | \&o) is shown as a function of the ramping time 
t r , along with the corresponding many-body overlap 
| (*(t r ) | *o) I- For the infinite Bose-Hubbard-model it 
has been shown [8( that t r 3> Vm ax /J 2 guarantees adi- 
abatic evolution. Our numerics of a finite-sized system 
agree with this result giving a near-unit overlap and small 
energy differences for ramping times of t r ~ 3Vb/ J 2 with 
Vo the initial interaction strength in the MI. 

Given that the ramping time t r is much longer than the 
time required to implement the MZI in Fig. [2b) we con- 
sider the influence of atom loss over a time t r . The build- 
up of long range correlations during the melting process 
has been found to be robust to atom loss processes 0. 
However, atom loss will cause uncertainty in the relative 
number of a and b atoms in the final state | $ t f ) ■ We use 



the model introduced in [Tl| to calculate the sensitivity 
scaling a assuming the state | Vt^f) is subject to atom loss 
at a rate k for time t r . As shown in Fig. |4tb) we find 
that a « 1 for nt r <C 0.1. In combination with the adi- 
abaticity condition on t r Heisenberg limited sensitivities 
can thus be achieved for n <C J 2 /30Vb. This condition 
is within the reach of current optical lattice technology 
where tunneling amplitudes J of a few hundred Hz and 
loss rates on the order of Hz can be achieved 0. Fur- 
thermore, additional fluctuations in the overall particle 
number N which occur from shot to shot do not affect 



of singly occupied lattice sites with atoms in the auxiliary 
state c. Additionally, an imperfect Raman transition dur- 
ing the creation of the initial MI state may leave a c atom 
pair within the | ^ a b) state. In both cases the presence 
of a small number of c-atom impurities will not destroy 
the coherence of the final SF state, as seen in recent ex- 
periments [2l| . Also these atoms will remain in state c 
throughout the whole process and thus not contribute to 
the particle number fluctuations 

In summary we have shown that quantum melting of a 
two-component MI in an optical lattice provides a viable 
route to engineer twin Fock on timescales t r smaller than 
current experimental atom-loss times. In particular the 
creation time t r does not scale with N. Our scheme also 
exploits the accurate controllability of two-atoms in a 
single lattice sites to minimize fluctuations in the relative 
numbers of a and b atoms which is crucial for achieving 
Heisenberg limited sensitivities cx 1/N. 
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the scaling a 17| . In the experiment a harmonic trapping 
potential is likely to be present and causes an outer shell 
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